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Abstract 

We report search results for the flavor-changing neutral current process b ^ d'y with a data 
sample containing 274 million B meson pairs accumulated at the T(45) resonance with the Belle 
detector at KEKB. We studied the exclusive decays B~ — p~7, -B° — p°7, and — uj, and find 
no significant signal. We set an upper limit for a combined branching fraction B{B — (p, 0^)7) < 
1.4 X 10~® at the 90% confidence level, which is normalized to that of B~ — assuming an 
isospin relation between the three modes. 

PACS numbers: ll.30.Hv, 13.40.Hq, 14.65.Fy, 14.40.Nd 
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The b ^ d'y process, shown in Fig. Cl^a), is a flavor changing neutral current transition 
that proceeds via loop diagrams in the Standard Model (SM). It is suppressed with respect 
to 6 — S7 by the Cabibbo-Kobayashi-Maskawa (CKM) factor \Vtd/Vts\'^ ~ 0.04 with a large 
uncertainty due to the lack of precise knowledge of Vtd- The exclusive modes B ^ pj 
and to'7, which are presumably dominant, have not yet been observed P, |^. They are 
also suppressed with respect to the corresponding exclusive decay B — > K*^ by |V^d/KsP, 
with corrections due to form factors, SU (3) breaking effects and the additional annihilation 
diagrams (Fig-TOb)), giving predicted branching fractions in the range (0.9-2.7) x 10~^ 
in the SM js], \^. Measurement of these exclusive branching fractions would improve the 
constraints on Vtd in the context of the SM, and would provide sensitivity to physics beyond 
the SM that is complementary to that from b —>■ 57. 

In this paper, we report the results of a search for the b ^ d'-f process using a data 
sample of (274 ± 3) million B meson pairs accumulated at the T(4S') resonance. The 
data areproduced in e~^e~ annihilation at the KEKB energy-asymmetric (3.5 on 8 GeV) 
collider and collected with the Belle detector j6| . The Belle detector is a large-solid-angle 
spectrometer that includes a silicon vertex detector (SVD), a central drift chamber (CDC), 
an array of aerogel threshold Cherenkov counters (ACC), time-of-fiight (TOF) scintillation 
counters, and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to identify muons (KLM). The dataset 
consists of two subsets with different inner detector configurations: for the first 152 million 
B meson pairs, a 2.0 cm radius beampipe and a 3-layer SVD were used; and for the remaining 
122 million B meson pairs, a 1.5 cm radius beampipe, a 4-layer SVD and a small-cell inner 
drift chamber were used |7|. 

We reconstruct the following final states: B^ —>■ p^7, p°7, and 5° u'y. (Charge 

conjugate modes are implied throughout this paper.) We also reconstruct control samples of 
B^ —>■ K*'^'~f and B^ —>■ K*^'~f decays. The following decay chains are used to reconstruct the 
intermediate states: tt^tt^, p° —>■ tt'^tt', uj tt+tt^tt^, K*~ — * K~Tr^, K*^ —>■ K'tt~^, 

and 7r° ^77. 

Photon candidates are reconstructed from isolated clusters in the ECL that have no 
corresponding charged track, and a shower shape that is consistent with that of a photon. 
The photon with the largest center-of-mass (CM) energy in the range 1.4 GeV < < 
3.4 GeV and in the barrel region of the ECL (33° < 9^ < 128° in the laboratory frame) 
is selected as the primary photon candidate. To suppress backgrounds from 7r° — 77 and 
— ;> 77 decays, we veto the event if the reconstructed mass of the primary photon and any 
other photon of 30 (200) MeV or more is within ±18 (32) MeV/c^ of the vr*^ (77) mass. These 
correspond to ±3(t windows, where a is the mass resolution. This set of criteria is referred 
to as the vr^/?] veto. For the primary photon, we sum the energy deposited in arrays of 3 x 3 
cells and 5x5 cells around the maximum energy ECL cell; if their ratio is less than 0.95, 
the event is vetoed. 

Charged pions and kaons are reconstructed as tracks in the CDC and SVD. Each track 
is required to have a momentum greater than 100 MeV/c and closest approach to the run- 
averaged interaction point within 2 cm in radius and ±5 cm along the z-axis, aligned op- 
posite the positron beam. We do not use the track to form the signal candidate if, when 
combined with any other track, it forms a candidate with a mass within ±10 MeV/c^ 
around the nominal mass and a displaced vertex that is consistent with a K^. We de- 
termine the pion (L^) and kaon [Lk] likelihoods from the ACC response, specific ionization 
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(dE/dx) in the CDC and TOF flight-time measurements for each track, and form a hke- 
hhood ratio C^/k = Lt,/{Lt, + Lk) to separate pions and kaons. We require Ct,/k > 0.85 
for pions, which gives an efficiency of 89% for pions and misidentification probabihty of 
~10% for kaons. For the ivy mode, we relax the requirement to C-jt/k > 0.8, which gives 
an efficiency of 94% for pions. (For kaons in the K*^ modes, we require C^/k < 0.4, which 
gives an efficiency of 76-80%.) In addition, we remove pion and kaon candidates if they are 
consistent with being electrons based on ECL, dE/dx and ACC information, or consistent 
with muons based on KLM information. 

Neutral pions are formed from two photons with invariant masses within ±10 (16) MeV/c^ 
of the mass, corresponding to a ~2cr (~3cr) window for p~7 and K*~'~f (ooy) modes. The 
photon momenta are then recalculated with a vr^ mass constraint. We require each photon 
energy to be greater than 30 MeV. We also require the CM momentum of the to be 
greater than 0.5 GeV/c for the p~7 and K*~'y modes. 

Invariant masses for the p and uj candidates are required to be within windows of 
±150 MeV/c^ (IF) and ±30 MeV/c^ (3.5F), respectively, around their nominal masses, 
where F is the natural width of each resonance. 

We form B candidates by combining a p or u candidate and the primary photon using 

two variables: the beam-energy constrained mass Mbc = \J {.E^ea.ml ^'^Y ~ IPb/cP and the 
energy difference Ai? = E*^ — E^^^^, where p*^ and E*^ are the measured CM momentum 
and energy, respectively, of the B candidate, and -Eboam the CM beam energy. The 
magnitude of the photon momentum is replaced by (-Ebeam ~ ^I/lo)!^ when the momentum 
p*Q is calculated. To optimize the event selection, we count Monte Carlo (MC) events in 
the region -0.10 GeV < < 0.08 GeV and 5.273 GeV/c^ < Mbc < 5.285 GeV/c^. We 
choose the selection criteria to maximize Ns / V^b, where Ns is the expected signal yield in 
this region assuming the branching fractions to be the SM value in Ref. j^, and Nb is the 
expected background yield in the same region. 

The dominant background arises from continuum events {e~^e~ — > qqij)) where the acci- 
dental combination of a p or candidate with a photon forms a B candidate. We suppress 
this background using an event shape discriminator JF, the B candidate polar angle 6'^ in 
the CM frame, the vertex separation Az, and the output of the i?-flavor tagging algorithm: 

1. The event shape discriminator JF j3is a Fisher discriminant constructed from 16 
modified Fox- Wolfram moments jsLllOj and the scalar sum of the transverse momentum 
of all charged tracks and photons. 

2. True B mesons follow a 1 — cos^6'^ distribution in cos^^, while candidates in the 
continuum background are uniformly distributed. 

3. In about 85% of events for the p°7 and wy modes, a fit can be successfully performed 
to determine the decay vertex of the candidate B meson as well as the origin of the 
remaining tracks in the event. The separation Az between these two vertices along 
the z-axis discriminates between continuum events, which have a common vertex, and 
signal events, whose decay vertices are displaced. 

4. The i?-flavor tagging algorithm described in Ref. returns the flavor of the other 
B meson (g = ±1), and a tagging quality r (0 < r < 1) which indicates the level 
of confldence in the flavor determination. The algorithm uses the particles in the 
event that are not associated with the signal B candidate, and provides additional 
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discrimination between signal, and continuum background where no true B meson is 
present. 

For each of the quantities cos^^ and Az, we construct one-dimensional hkelihood 
distributions for signal and continuum. Signal distributions are modeled with an asymmetric 
Gaussian function for JF, |(ao — a2COS^^g) for cos^|j, and an exponential convolved with a 
Gaussian resolution function for Az; continuum background distributions are modeled with 
an asymmetric Gaussian function for {bo — b2COs'^9%) for cos^^, and a sum of three 
Gaussian functions with a common mean for Az; the coefficients Oq, 02, ^'o are close to unity 
while the coefficient 62 is close to zero. 

Since JF, cos^^ and Az are independent quantities, we form a likelihood ratio TZ = 
Cs/{Cs + C-c) to combine them, where Cg and Cc are products of cos6'|j (and Az if 
available) likehhood distributions for signal and continuum, respectively. The likehhood 
distribution for the background Az distribution is determined from data in the sideband 
region 5.20 GeV/c^ < Mbc < 5.24 GcV/c^, |A£;| < 0.3 GeV; aU the other hkelihood distri- 
butions are determined from MC samples. 

We determine the Az likelihood function separately for the two datasets as the vertex res- 
olution is improved in the SVD2 with respect to the SVDl. As a consequence, we introduce 
two sets of likelihood ratios and selection criteria for each decay mode. 

On the plane defined by the tagging quality and likelihood ratio, (r, 7^), signal tends to 
populate the edges at r = 1 and TZ — 1] continuum tends to populate the edges at r = 
and TZ — Q. We select the events in a signal enriched region defined hy TZ> TZi for r > ri, 
and TZ > 1 — «(1 + r) for r2 < r < ri, where the parameters ri, r2, -Ri and a are mode 
dependent and are determined so that Ns / \/ Ns + Nb is maximized (we use this quantity 
instead of Ns / \J Nb because of the limited statistics of the MC simulation sample that is 
used in this procedure). The values are ri — 0.85, r2 = 0.01, Ri > 0.8 and a — 0.025 for 
the ^"7 mode in the first subset of data; similar values are used for the other subset and for 
the p~7 and wy modes. We define the rest of the the background enriched region. 

We consider the following B decay backgrounds: B — > K*j, other B — > Xgj processes, 
B — > p7r° and a;7r°, B pr] and ur], B~ P~P^, other charmless B decays, and b ^ c 
backgrounds. We find the b ^ c background to be negligible. The B — > K*'y background 
may mimic the signal decay B ^ p'-f if the kaon from K* is misidentified as a pion. In order 
to further suppress B —>■ K*'y, we calculate Mktv, where the kaon mass is assigned to one 
of the pion candidates, and reject the candidate if M^tt < 0.96 (0.92) GeV/c^ for the p^^ 
{p--f) mode. The decay chain W K*^-t, isT^Tr", tt+tt" has the same final 

state as B^ — > 00^, and has a small contribution due to the tail of the K* Breit-Wigner line 
shape. In addition, B — * K*'~f and other B Xsj decays contribute to the background 
when the p and u candidates are selected from a random combination of particles. Charmless 
decays with a tt*^ or 77 ^ 77, B pir^, cutt^, prj and urrj, may mimic the signal if one of the 
photons from 7r° or rj decay is soft and undetected by the vr^/ry veto condition. To suppress 
this background, we calculate the cosine of the helicity angle ^hci, and reject the candidate 
if I cos 6*1161 1 > 0.8 (0.6) for the p°7 and luj {p~j) modes. Here, 6*1161 is defined as the angle 
between the tt"*" and B momentum vectors in the p rest frame, or the angle between the 
normal to the u decay plane and the B momentum vector in the cu rest frame. The decay 
B~ P~P^, p~ 7r~7r° also contributes to the B^ — > p°7 mode when both one pion 
from the p~ decay and one photon from the 7r° decay are soft and undetected. The other 
charmless decays have small contributions and are considered as an additional background 
component when we extract the signal yield. 
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The reconstruction efficiency for each mode is defined by the fraction of the signal yield 
remaining after all selection criteria, where the signal yield is determined from an extended 
unbinned maximum likelihood fit to the MC sample. The signal distributions are mod- 
eled as the product of a Gaussian function for Mbc and an empirical function known as 
the Crystal Ball line shape 0] to reproduce the asymmetric ECL energy response for AE". 
The background c omp onent is modeled as the product of a linear function for AE, and an 
ARGUS function 'IS] for Mbc- From the fit, we find efficiencies of about 5% as listed in 
Table m The systematic error on the efficiency is the quadratic sum of the following contri- 
butions estimated from control samples: the photon detection efficiency (2.2%), measured 
from radiative Bhabha events; charged tracking efficiency (1.0% per track) from partially 
reconstructed D*^ — > Z)°7r+, —>■ fC^vr+vr", Kg — »• 7r"'"(7r~); charged pion identification 
(1.0% per pion) from D*^ —>■ D^n^, K^n^; neutral pion detection (4.6-7.3%) from rj 

decays to 77, h'^tt^tt^ and Stt*^; TZ-r and tt^/t] veto requirements (5.4%) from B'^ —>■ D^n^, 
K-n+] and MC statistics (0.9-1.5%). 

We perform an unbinned maximum likelihood fit to the data in the (Mbc AE) fit region 
bounded by \AE\ < 0.3 GeV and Mbc > 5.2 GeV/c^, simultaneously for the three signal 
modes (collectively referred to as i? — (p, uj)'y) and the two B — >• K*j modes. We fit the 
two data subsets simultaneously, so that in total ten distributions are included in the fit. 
We define the combined branching fraction B{B {p,uj)'~f) = B{B~ P~l), assuming the 
isospin relation Q B{B- p-7) = 2^B{W p°-f) = 2^B{W ^ uj-f), where we use 

^ = 1.086 ± 0.017 0. We also assume B{B K*-f) = BiB' K*--f) = ^B{W 

#%). 

We describe the events in the fit region using a sum of the signal, continuum, K*^, and 
other background hypotheses. Parameters of the signal description for Mbc are calibrated 
using B^ D^Ti^ samples. Those for AE are calibrated using the result of a fit to the 
AE distribution of the B K*'-/ control sample, in which the mean and the width are 
allowed to fioat. We use the branching fractions for B — > (p, u;)7 and B K*^ as the 
free parameters, from which the signal yield for each channel is deduced using the efficiency 
for each mode. We assume the efficiency systematic errors are fully correlated when we 
evaluate the systematic error. The continuum background is modeled as the product of a 
linear function for AE whose slope is allowed to fioat, and an ARGUS function for Mbc whose 
parameters are fixed from a comparison between data and MC in the background enriched 
region and MC in the signal enriched region. The continuum contribution in the data sample 
is allowed to fioat. The size of the K*'y background component in each {p,uj)'y channel is 
constrained using the fit to the K*'y events. The contributions of the other backgrounds are 
fixed using known branching fractions or upper limits. The free parameters in the fit are 
therefore the branching fractions for B — >■ (p, a;)7 and B K*'y, five continuum fractions, 
and five continuum AE slopes. 

We also perform individual fits to the three signal modes and the two B K*'j modes. 
The two data subsets are fitted simultaneously for each mode. The size of the K*^^ back- 
ground in the p°7 mode and K*~'-f background in the p~7 mode are fixed according to the fit 
results to the B —>■ K*^ modes and the known particle misidentification probabilities. The 
free parameters in each fit are the signal yield, the continuum fraction, and the continuum 
AE slope. 

Results of the simultaneous fit are shown in Fig.|21and given in Table|l] The simultaneous 
fit gives a significance of 1.9 standard deviations, where the significance is calculated as 
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y'— 2 ln(£o / '^max) , and £max (-^^o) is the maximum likelihood from the fit when the signal 
branching fraction is floated (constrained to be zero). In order to include the effect of possible 
systematic error in the significance calculation, we change each parameter by one standard 
deviation in the direction that gives the smallest resulting significance. The systematic 
error in the signal yield is estimated by varying each of the fixed parameters by its standard 
deviation, and then taking the quadratic sum of the deviations in the signal yield from the 
nominal value. The combined branching fraction is B{B (p,uj)'j) = {0-72^q '^1~^q'^j) x 
10~^, where the first and second errors are statistical and systematic, respectively. 

Since the significance is small, we quote a 90% confidence level upper limit Bgo using 
the relation Jq^" C{x)dx = 0.9 C{x)dx, where C{x) is the likelihood function with the 
branching fraction fixed at x. The systematic error is taken into account assuming a Gaussian 
distribution. We find 

B{B ^ {p,uj)-f) <lAx 10-^ (1) 

at the 90% confidence level. Results of the fits to the individual modes are also given in 
Table HI 

A similar fit procedure is performed by using the ratio of branching fractions B{B — »• 
{p,uj)j)/B{B K*'y) instead of B{B {p,uj)'y), so that the systematic error partially 
cancels. We find B{B {p,uj)-f)/B{B K*-^) < 0.035 at the 90% confidence level. One 
can use this result to constrain Vtd- for example, using the prescription given in Ref. [l5l |. 

g(i?^(p,^)7) 
B{B K*j) 

where the form factor ratio ( = 0.85 ± 0.10 and SU{3) breaking effect AR = 0.1 ± 0.1, 
we obtain |V((i/Vjs| < 0.21 at the 90% confidence level. This limit is consistent with other 
determinations of |Vtd/Vts| jl^ . 

In conclusion, we search for the b ^ d'j process using a simultaneous fit to the B — >■ p7, 
B^ — > u'j and B K*'y modes. The upper limit we obtain is already within the range of 
SM predictions 0, |^ and can be used to constrain Vtd- 
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(a) loop diagram (b) annihilation diagram 




FIG. 1: (a) Loop diagram for b ^ d'j and (b) annihilation diagram only for B — > p 7. 
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FIG. 2: Projections of the simultaneous fit results to Mbc (in the region —0.10 GeV < AE < 
0.08 GeV) and AE (in the region 5.273 GeV/c^ < Mbc < 5.285 GeV/c^) for the individual modes 
and the sum of them. Lines represent the total fit result (solid), signal (dashed), continuum 
(dotted), and B decay background (dot-dashed) components. 



TABLE I: Efficiencies, significances, and 90% confidence level upper limits for the branching 
fractions. 



Mode 


Efficiency 


Significance 


Upper limit 




(zbsyst) 




(90% C.L.) 


B (p, combined 




L9 


1.4 X 10"^ 


B- p-j 


(5.5 ± 0.4)% 


2.1 


2.2 X 10"^ 


W ^ 


(3.9 ± 0.3)% 


0.6 


0.8 X 10-6 




(3.9 ± 0.4)% 


0.2 


0.8 X 10-6 
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